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ABSTRACT 

New photometry of RRab and RRc stars in 10 Cen is used to calibrate their 
absolute magnitudes My as a function of: a) metallicity; and b) the Fourier parameters 
of light curves in the V band. The zero point of both calibrations relies on the distance 
modulus to the cluster derived earlier by the CASE project based on observations of 
the detached eclipsing binary OGLE GC17. For RRab variables we obtained a relation 
of My = (0.26 ± 0.08)[Fe/H] + (0.91 ± 0.13). A dereddened distance modulus to the 
LMC based on that formula is /j,q — 18.56 ± 0.14 mag. The second calibration of My, 
which is based on Fourier coefficients of decomposed light curves, results in the LMC 
distance of /iq = 18.51 ± 0.07 mag. 
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1 INTRODUCTION 

, The Cluster AgeS Experiment (CASE) is a long-term 
project whose main goal is to determine the distances and 
ages of nearby globular clusters (Thompson et al. 2001). The 
primary distance indicators used by the CASE group are de- 
tached eclipsing binaries. The search for these relatively rare 
objects is performed using 1-meter class telescopes. Usually 
for each target cluster, several hundred CCD frames are col- 
lected over one or more observing seasons. Suitable candi- 
dates are then observed with larger telescopes (eg. Thomp- 
son et al. 2001). A secondary benefit of the project's survey 
phase is the detection of large samples of variables of various 
types including RR Lyr and SX Phe stars. 

The first globular cluster for which the CASE project 
determined a distance was lo Centauri. An analysis of pho- 
tometric and spectroscopic data for an eclipsing binary 
OGLE GC17 yielded an apparent distance modulus equal 
to (m - M) v = 14.09 ± 0.04 mag (Kaluzny et al. 2002). 

Time series photometry of lo Cen obtained by the CASE 
project in the 1999 and 2000 seasons resulted in light curves 
of almost 400 variables from the cluster field (Kaluzny et al. 
2003). The most numerous group were the RR Lyr stars. The 
new precise and well-sampled BV light curves of these stars 
can be supplemented with information about their metallic- 
ities (Rey et al. 2000). In this paper we use these data to 
calibrate different formulae which permit the calculation of 
absolute magnitudes of RR Lyr variables. 



2 M v - [FE/H] RELATION 

It has been known for many years that the absolute magni- 
tude of RR Lyr stars My is a function of metallicity [Fe/H] 
(Sandage 1981a,b). This function is often assumed to be a 
simple linear relation of the form: 

M v = a[Fe/H] + /3 (1) 

There is a lack of general consensus on the exact values 
of the a and (5 parameters. The most extreme values of the 
slope of the relation (1) were given by Sandage (1981a,b) 
with a = 0.35 and Fusi Pecci et al. (1996) with a = 0.13 ± 
0.07. The estimates of the zero point of the relation vary 
from around 1.1 mag (Gould & Popowski 1998) to 0.8 mag 
(Gratton et al. 1997, Caloi et al. 1997). 

It is worth pointing out that some theoretical models 
predict that a simple linear relationship of form (1) may not 
exist at all. This is caused by the fact that even variables 
of the same metallicity have different luminosities depending 
on the direction of their evolution across the instability strip 
(Lee 1991). More complex relations between My of RR Lyr 
stars and their metallicities were documented, for example, 
by Caputo (1997). 

There is also observational evidence for a non-linear 
dependence of My on [Fe/H]. For example, Caputo et al. 
(2000), analyzing RR Lyr variables from Galactic globu- 
lar clusters, obtained a = 0.17 ± 0.04 for [Fe/H] < -1.5 
and a = 0.27 ± 0.06 for [Fe/H] > -1.5. The systematic er- 
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rors caused by using the simple linear relation (1) may be 
minimized by excluding highly evolved RR Lyraes from the 
analyzed samples. In particular, the formula should be used 
with caution when analyzing variables belonging to globular 
clusters showing "blue" horizontal branches. 

Recent investigation of the Mv — [Fe/H] relation in ui 
Cen was done by Rey et al. (2000). They used [Fe/H] metal- 
licities derived from the hk index of the Caby photometric 
system. Their intensity averaged magnitudes of RR Lyr stars 
were taken from photographic photometry from Butler et al. 
(1978) and the CCD photometry of Kaluzny et al. (1997). 

The new V-band light curves of RR Lyr variables in ui 
Cen obtained by the CASE (Kaluzny et al. 2003) have about 
three times as many observed points as the photometry of 
Kaluzny et al (1997). This permits a more precise estimation 
of mean < V > magnitudes and allows the elimination of 
objects with unstable light curves from the calibration sam- 
ple. Adopting an apparent distance modulus of the cluster 
as determined by (Kaluzny et al. 2002) we can obtain a new 
Mv — [Fe/H] relation for RR Lyr stars in cu Cen. 

The sample used for the calibration includes 122 stars 
with stable light curves of good quality and metallicities de- 
termined by Rey et al. (2000). We decided to remove the 
RRab variable V52 from the further analysis. Its high lumi- 
nosity suggests that it is a foreground star. Additionally, van 
Leeuwen et al. (2000) gives only a 45% membership proba- 
bility for this variable, based on a proper motion study. 

The Mv — [Fe/H] dependence for remaining 121 vari- 
ables is shown in the upper panel of Fig. 1. The middle panel 
shows the same relation but for a sub-sample consisting of 
67 RRab variables. The errors in [Fe/H] were taken directly 
from Table 5 of Rey et al. (2000) . For objects with unknown 
error of [Fe/H] we assumed it to be equal to 0.2 dex as sug- 
gested by Rey et al. (2000). 

The intensity averaged magnitudes < V > obtained 
from light curves containing 500-800 points have internal 
errors smaller than 0.001 mag. We estimate that the external 
error of the zero point of the photometry is about 0.02 mag. 

Transformation from < V > to Mv is performed using 
the distance modulus whose error is 0.04 mag. Thus, combin- 
ing this value with 0.02 mag error resulting from computing 
the mean magnitudes < V >, we assumed that individual 
absolute magnitudes of RR Lyr stars from our sample have 
uncertainties of 0.05 mag. 

Knowing the errors we were able to fit straight lines to 
the graphs shown in Fig. 1. For the sample of 121 RR Lyr 
stars we obtain the following relation: 



-.5 



M v = (0.21 ± 0.05) • [Fe/H] + (0.77 ± 0.08) 
and for 67 RRab variables is in the form: 
M v = (0.28 ± 0.07) • [Fe/H] + (0.89 ± 0.12) 



(2) 



(3) 



The values of a derived above may be biased by the 
presence of extremely evolved objects in our sample. The 
evolutionary models of Lee (1990) suggest that the RR Lyrac 
stars in clusters having a very blue HB and with a metal- 
licity in the range of —2.0 < [Fe/H] < —1.6 are highly 
evolved stars. They have significantly brighter magnitudes 
and longer periods than those near the zero age horizontal 
branch (ZAHB). 

Rey et al. (2000) investigated this problem in uj Cen. 
Analyzing the period-amplitude relations for different ranges 
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Figure 1. The relation between My and [Fc/H] for: a) the whole 
sample; b) sub-sample of RRaft variables; and c) for unevolved 
RRafe variables. The solid lines arc best linear fits. 



of metallicities they discovered that a significant sample of 
evolved RR Lyr stars exists only for variables with —1.9 < 
[Fe/H] < —1.5. Looking at their Fig. 9b one can clearly see 
that these evolved objects are also characterized by periods 
longer than 0.7 days. 

Thus we decided to exclude from our sample all RRafe 
variables having metallicities in a range of —1.9 < [Fe/H] < 
— 1.5 and pulsation periods longer than 0.7 days. Our final 
sample consists of 53 RRab variables shown in the lower 
panel of Fig. 1. The resulting relation for Mv is as follows: 



M v = (0.26 ± 0.08) ■ [Fe/H] + (0.91 ± 0.13) 



(4) 



The slope of this relation falls inside the 0.13-0.35 range 
presented in the literature. Recently, Chaboyer (1999) in his 
review of globular cluster distance determinations adopted 
a = 0.23 ± 0.04 arguing that 1-a range of this value encom- 
passes the majority of recent determinations for a. 

Our results do not differ significantly from those pre- 
sented by Rey et al. (2000), due to the good agreement of the 
zero points of the CASE photometry, as well as the CCD and 
photographic photometry used by Rey et. al. (2000). The 
mean difference in V magnitude between both data sets con- 
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sisted of 53 RRa6 stars is < Vcase- V Rey >= 0.007±0.010 
mag. 

Having established the relation (4), we can use it for 
determining the distance modulus to the Large Magellanic 
Cloud - (ilmc- We used the data of Walker (1992), who 
summarized the CCD photometry of 182 RR Lyr variables 
belonging to seven LMC globular clusters. We decided to 
exclude from our analysis variables from the cluster NGC 
1841 because it is most probably located about 0.3 mag 
closer than the main body of the LMC. The mean value of 
the reddening free magnitude < Vb > of 160 RR Lyr stars 
from the remaining six clusters is 18.98 ± 0.03 mag. The 
mean metallicity of this sample is [Fe/H] = —1.9 ± 0.2. 

Knowing the observed magnitude of RR Lyr stars from 
the LMC and computing their absolute magnitude using 
relation (4) we determined that the distance modulus to 
the LMC is equal to [ilmc = 18.57 ± 0.20 mag. This result 
supports the "long" distance scale. 

Most recently, Clementini et al. (2003) presented new 
photometry and spectroscopy for more than a hundred RR 
Lyrae stars in two fields located close to the bar of the LMC. 
The dereddened apparent average luminosity of the vari- 
ables from their sample was < Vb >= 19.064 ± 0.064 mag 
at the mean metal abundance of [Fe/H] = —1.48 ± 0.03. 
Their metallicity was tied to the Harris (1996) scale which 
differs by about 0.06 dex from the Zinn & West (1984) scale 
used by Rey et al. (2000) . Thus the mean metallicity of their 
sample is [Fe/H] = -1.54 ±0.03 on the Zinn & West (1984) 
scale. 

Application of equation (4) for their data gives the dis- 
tance modulus to the LMC equal to [ilmc = 18.55 ± 0.19 
mag. 

Combining the results of photometry obtained by both 
Walker (1992) and Clementini et al. (2003) we find that our 
My-[Fe/H] relation gives ^lmc = 18.56 ± 0.14 mag. 



3 ABSOLUTE MAGNITUDE OF RR LYR 

STARS AS A FUNCTION OF THE FOURIER 
PARAMETERS 

3.1 RRat stars 

3.1.1 Methods 

On the one hand, nonlinear pulsation models of RR Lyr stars 
suggested that the luminosity of these stars may be uniquely 
related to the pulsation period and the shape of their light 
curves (e.g. Simon & Clement 1993). This strongly suggests 
that there is a close correlation between the Fourier coeffi- 
cients of the light curves and their corresponding periods and 
luminosities. Such correlations could potentially improve the 
distances derived for these stars. 

On the other hand serious attempts to find such em- 
pirical relations pioneered by Kovacs & Jurcsik (1996, here- 
after KJ96) and most recently updated by Kovacs & Walker 
(2001) - hereafter KW01 - bore mixed success. The empirical 
procedure is essentially reduced to adopting a linear combi- 
nation of the period and M selected Fourier parameters T 
for predicting the RR Lyr absolute magnitude: 

M 

M V {PF) = C + CP log P + C^Fm (5) 

m=l 



The coefficients c m and cluster distance moduli dk are fitted 
to minimize scatter \ 2 m the apparent average magnitudes 
Vkn summed over N k RR Lyr stars in K clusters: 

K N k 

X = J2J2[Vkn -d k - M V (PF)] 2 (6) 

fc = l n=l 

Problems are indicated by an apparent lack of convergence 
from the whole procedure to a unique solution. Different 
data sets yield substantially different relations, both in re- 
spect to the optimum parameter set and in the numerical 
values of the involved coefficients. For RRab stars KJ96 ob- 
tained: 

M v = 1.221 - 1.396P - 0.477Ai + 0.103<^ 3 i (7) 

while KW01 using 383 RRafe stars from 20 globular clusters 
recommend either of three: 

M v = -1.820 log P- 0.805.4 1 +c (8) 
M v = -1.876 log P - 1.158Ai + 0.82L4 3 +c (9) 
M v = -1.963 log P- 1.124Ai + 0.830A 3 + 

+ 0.011051 +c (10) 

Note that the sine phase convention is used here unless 
otherwise stated. A worrying inconsistency is shown, for ex- 
ample, by the values of cp and cai in Eq. (7) and (8). 

Several different causes might be responsible for these 
problems: 

(i) The original problem connected to the sensitivity of 
high order Fourier coefficients to poor data melted away with 
abundance of the high quality observations (e.g. Kaluzny et 
al. 1997, KW01 and the present work). 

(ii) The more persistent problem is connected to the ap- 
parent inability of the procedure to account for the hidden 
diversity of properties related to metallicity, evolutionary 
status and possibly misidentified pulsation mode, all capa- 
ble of producing several different trends in the data. This 
problem is acknowledged by rejecting outlying stars from 
the rest by the D m distance criterion (see e.g. Kovacs & 
Kanbur 1998). The risk is that the omitted stars hide much 
information on the nature of the physics involved. 

(iii) One reason for widely different values of the coef- 
ficients obtained from the different data sets could be the 
strong correlation between the Fourier coefficients. Such fits 
could be sensitive to a group of coefficients while changing 
loads between individual members of the group matters lit- 
tle. 

(iv) Finally, some evidence discussed further in this paper 
indicates that magnitudes of RRab stars might also depend 
on a factor, so far unaccounted for, independent of the period 
and shape of the light curve. We recommend future trials 
with colors sensitive to temperature and/or metallicity. 

In order to clarify the situation and possibly to identify 
the cause of these problems we undertook the same kind of 
analysis from scratch and for an entirely new set of data. 
Our present sample consists of 76 RRafo variables observed 
in u> Cen by Kaluzny at el. (1997) and Kaluzny et al. (2003). 
The photometry in these two papers is of very good qual- 
ity and shows no systematic differences in the zero point 
(see Kaluzny et al. 2003 for details). Following recommen- 
dations of Kovacs & Kanbur (1998) for analysis we used only 
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N^cen = 56 stars with D m < 3. A clear advantage of our 
data is very consistent photometry obtained within the same 
project and using the same telescope. The vast majority of 
stars was observed at least 500 times, hence we obtained 
very precise and reliable values of the Fourier coefficients. 

To minimize errors in the coefficients and their 
correlation we converted our projection into orthogo- 
nal trigonometric polynomials (Schwarzenberg-Czerny f 996, 
Schwarzenberg-Czerny & Kaluzny f996). As all stars come 
from the same cluster we do not suffer from the distance in- 
determinacy and differential interstellar reddening because 
there are no systematic differences in E(B — V) in the field 
of uj Cen (Schlegel et al. f998). Additionally, because of the 
significant spread in metallicities within u> Cen, our sam- 
ple is excellent for studying metallicity related intra-cluster 
differences unaffected by distance differences. 

Figure 2 shows the amplitudes Aj, amplitude ratios Rji 
and phase combinations <pji as the function of the period P 
for all RRab variables from our sample. Comparing them 
with the same relations for the KWOf sample (see their Fig. 
f) we conclude that our relations are significantly tighter 
and correlate very well with the period. 

The interesting thing we noted is the clear change of 
the slope in the <j>si - P and </>6i - P relations occurring 
at P w 0.75 days. Arranging the RRab stars from to Cen 
according to the increasing period one can see that it is 
connected with a disappearance of the "bump", clearly seen 
at phase around 0.7 in the light curves of variables with 
periods shorter than 0.75 days. 

The basic properties of our sample of KRab stars from 
lu Cen are collected in Table f . 

3.1.2 Results for u) Cen 

We performed several fits involving P and up to 3 different 
Fourier phase and amplitude coefficients drawn from the first 
5 harmonics. To remove the constant term indeterminacy, we 
adopted as fixed the distance modulus 14.09±0.04 
mag of Kaluzny et al. (2002). The results are listed in Ta- 
ble 2. We explicitly list the fitted formula and its standard 
deviation T>. To test the numerical self consistency of our 
fits we use each formula to recalculate /i^cen and its stan- 
dard deviation a. These are also listed in Table 2. Its hardly 
surprising that V/a w s/N^Cen- We shall refer to different 
fitted formulae by their consecutive numbers in the Table 2, 
from (Fl) to (Fll). 

Inspection of Table 2 shows that no formulae involving 
only amplitudes (e.g. Fl & F5) are satisfactory. This is in 
marked contradiction to expectations from Eqs. (8) and (9). 
It is particularly surprising as the latter formula performed 
best on the data from KWOf. We observe a strong correla- 
tion effect of type (iii) in Sect. 3.1.1. In Fl & F5 only the 
amplitude coefficients vary while the rest remains remark- 
ably stable, and the quality of the fit (D) is unaffected. In F5, 
contributions from Ai and A3 simply cancel out to the value 
in Fl for A\ alone. Most other fits involving both amplitude 
and phase appear satisfactory. It seems that the set of pa- 
rameters originally selected by KJ96 performs quite well in 
our case, except that our coefficients in F3 are completely 
different from those in Eq. (7). Formally, our best fit is F6. 
However, the improvement of \ 2 f° r F6 compared to that 
for F3 yields a Fisher-Snedecor test value of F(l, 52) = 7.85, 



only marginally significant at the level of 99%. At the face of 
these values one could say that F2, F3, F6-F10 yield fits of 
similar quality and inclusion of more parameters are hardly 
justifiable. It looks like the whole procedure is incapable of 
yielding better accuracy than these latter formulae. 

In order to reveal any remaining trend in the residu- 
als from the fit in each panel of Fig. 3 we plot for all stars 
their absolute magnitude My predicted by a given formula 
against average apparent magnitude < V >. Dots corre- 
spond to the formulae from Table 2 and circles to Eqs. (8- 
10) from KW01, the latter shifted by a constant. The lines 
mark ideal relations My =< V > —fi^cen- Inspection of 
Fig. 3 reveals that compared to Eqs. (8-10) our formulae 
involving phase perform particularly well for brighter stars. 
The common feature in all panels is the markedly lower in- 
clination of the trend in the points than the ideal relation. 
This is worrisome as apparently the fitted formulae are un- 
able to reproduce the observed span of magnitudes. This 
effect is also present in similar attempts to predict magni- 
tudes of cepheids using shapes of their light curves (Ligeza 
& Schwarzenberg-Czerny 2000). In our opinion the presence 
of such an effect might indicate that some hidden factor 
(different from the period and shape of the light curve) sig- 
nificantly influences the average magnitudes of RRab stars 
(see point iv in Sect. 3.1.1). 

3.1.3 Application to LMC 

To verify the usefulness of the formula F6 we decided to de- 
termine the distance modulus of the LMC. To do this we 
required high quality photometry of RR Lyr variables from 
the LMC. The first source we verified was the photometry 
of RR Lyr variables in seven globular clusters placed in the 
LMC (Walker 1992). Unfortunately, these light curves con- 
tained only about 30 points, which is insufficient to obtain 
valuable information about high order Fourier coefficients. 

The photometry of 68 000 variable stars in the LMC 
was performed by the OGLE team and recently published 
by Zebrun et al. (2001). Unfortunately, due to the fact that 
OGLE searches for microlensing phenomena, the vast ma- 
jority of their data is collected in the /-band. Additionally, 
the exposure times are too short to produce very high qual- 
ity data on RR Lyr variables. The photometry described by 
Zebrun et al. (2001) was collected in 21 fields and for most 
of them the number of U-band frames was around 30-40. 
Only four fields, namely SC2, SC3, SC4 and SC5, were ob- 
served over 50 times in V and thus the light curves of RR 
Lyr stars located in these fields have photometry suitable 
for our tests. 

The total number of RRafe stars in the above-mentioned 
four fields of the LMC was 1428. Only 13 of them have 
light curves characterized by the parameter D m < 3 and we 
used our F6 equation only for these. Such a large number of 
rejected stars cannot remain without discussion because it 
could suggest that our formulae are incapable of predicting 
the behavior of LMC RRab stars. The OGLE observations of 
RR Lyr variables are only a by-product of a project designed 
for other purposes, thus the number of V filter observations, 
their exposure length and distribution are not optimal for 
measuring of such fine effects as the high harmonics and 
amplitudes of the light curves of RRab stars. 

Leaving aside the problem of a large number of rejected 
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Figure 2. The amplitudes Aj, amplitudes ratios Rji and phase combinations ipji in the function of the period P for all RRafe variables 
from our sample. 



stars, we proceeded to apply the remaining 13 stars to deter- 
mine the LMC distance modulus. As we are concerned with 
the zero point only, the problems with the My modeling 
and intrinsic scatter of stars pose less danger. The resulting 
My — < V > relation for these 13 variables is shown in Fig. 
4. Again the solid line has a slope of unity. Averaging the 
differences between the observed and absolute magnitudes 
of the stars from this sample we obtain the value of the ap- 
parent distance modulus to the LMC as equal to 18.83±0.04 



mag, where the error is a standard deviation of the mean of 
13 estimates. 

In order to obtain the true distance modulus to the main 
body of the LMC, knowledge about interstellar reddening 
is needed. Due to the significant angular size of the LMC 
in the sky, the reddening varies depending on the position. 
However, the OGLE fields which we use are near the center 
of the LMC and adopting the average value of E(B — V) 
taken from recent estimates of reddening should not intro- 
duce large systematic errors. As was shown, for example, by 
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Figure 3. The My — < V > relations for the different formulae connecting the absolute magnitudes of RRab stars with their light curve 
parameters. Solid lines in all panels have a slope of unity. Filled circles correspond to My values computed from our formulae and open 
circles to the formulae of KW01. The My values obtained from KW01 equations are shifted by 0.4 mag for clarity. 
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Table 1. Basic parameters of RRab in ui Cen based on the photometry of Kaluzny ct al. (1997) and Kaluzny et al. (2003). 
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Alves et al. (2002), Fitzpatrick et al. (2002), Dutra et al. 
(2001) and Groenewegen & Salaris (2001) the interstellar 
reddening toward the LMC varies from 0.086 to 0.12 mag. 
Thus we simply assume that E(B — V) is equal to 0.10±0.02 
mag and our estimate of the true distance modulus to the 
LMC is hlmc = 18.51 mag. Taking into account the stan- 
dard deviation of the mean value of our 13 estimates, the 



error in the absorption in V-band introduced by the error in 
E(B — V) and the error in calibrating the zero point from 
the relation F6 of our final estimate of the distance to the 
LMC is fiLA4C = 18.51 ± 0.07 r ± 0.04 s mag. 
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Figure 5. The amplitudes Aj, amplitudes ratios Rji and phase combinations <f>ji in the function of the period P for all RRc variables 
from our sample. 
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Table 2. Formulae for My of RRafc stars with their T> parameters and the resulting distance moduli to u> Ccn with its standard deviations 



No 


Equation 










V 




a 


i 


ivi y — 


fl 979 


l.oDo lOg r — 








n noon 
U.Uyzu 


1 a noon 


U.ulzo 


2 


My — 


— U.y ( U 


- 3.055 logP 


+ 0.734Ai 


+ 0.333041 




0.0692 


14.0903 


0.0092 


3 


My = 


-2.858 


- 3.016 log P 


+ 0.529Ai 


+ 0.488</>3i 




0.0722 


14.0917 


0.0096 


4 


My = 


-0.911 


— z.ovi- log i 


i n Q07 A 
+ V.oZ l /±i 


+ O.163051 




U.D/oo 




0.U1U4 


5 


My = 


0.286 - 


1.872 logP - 


0.717Ai + 


O.4IOA3 




0.0928 


14.0895 


0.0123 


6 


M V = 


-1.000 


- 3.155 logP 


+ 0.095Ai 


+ 1.698A 3 - 


f 0.356041 


0.0673 


14.0907 


0.0088 


7 


My = 


-0.822 


- 3.114 log P 


+ 0.756Ai 


+ 0.453041 


- 0.082051 


0.0690 


14.0899 


0.0090 


8 


My = 


-3.135 


- 3.162 logP 


- 0.211Ai 


+ 2.OO2A3 - 


f 0.542031 


0.0696 


14.0902 


0.0091 


9 


M v = 


-1.214 


- 3.068 log P 


+ 0.729/1! 


+ 0.058031 


+ 0.298041 


0.0699 


14.0918 


0.0092 


10 


My = 


-2.873 


- 3.017 log P 


+ 0.528Ai 


+ 0.493031 


- 0.002031 


0.0730 


14.0899 


0.0096 


11 


My = 


-0.946 


- 2.666 log P 


- 0.219Ai 


+ 1.413A 3 + O.17405i 


0.0776 


14.0916 


0.0102 



M, 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



<V> 

Figure 4. The My- < V > relation for the 13 OGLE KRab 
variables from the LMC. The solid line has a slope of unity. 



3.2 RRc stars: near sinusoidal pulsators 

3.2.1 Is the lurninosity-shape correlation real for RRc? 

In fact theoretical investigation of correlations among phys- 
ical parameters of stars and their pulsation light curve for 
RRc stars started before those for RRab. Investigations by 
Simon & Teays (1982) and Simon (1989) relying on a large 
number of hydrodynamic pulsation models indicated that 
L correlates with (^31 and predominantly with P. Simon & 
Clement (1993) came out with the following formula for RRc 
luminosity: 



log L — 1.04 log P - O.O580|i + 2.41 



(11) 



where 03i indicates cosine phase convention, differing from 
our sine convention by it. The form of this formula resembles 
Eq. (7) except for a linear transformation of units. 

Our sample contains 54 RRc stars from w Cen observed 
by Kaluzny et al. (1997, 2003) within the OGLE and CASE 
projects. For this sample we selected only stars free of any 
complications, such as period changes, multiple periods, or 
low-amplitude noisy light curves. For the RRc stars we fol- 
lowed the same procedure as for the RRab. First we de- 
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Figure 6. The My— < V > relations for the different formulae 
connecting the absolute magnitudes of RRc stars with their light 
curve parameters. Solid lines in all panels have a slope of unity. 
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termined their Fourier parameters and then proceeded to 
fit formulae akin to Eq. (5) by minimizing \ 2 ■ On the one 
hand, the quality of our photometry for RRc is as good as 
ever. On the other hand, the small amplitude and nearly si- 
nusoidal light curve of RRc stars produced increased errors 
for high harmonics, so that above the 6-th harmonic they 
exceeded 10% of the value. The periods, magnitudes and 
Fourier parameters of the RRc stars from our sample are 
listed in Table 3. The amplitudes Aj, amplitudes ratios Rji 
and phase combinations <f>ji in the function of the period P 
for all RRc variables from our sample are shown in Fig. 5. 

These data were fitted with formulae of the type of Eq. 
(7), using P and up to 3 Fourier coefficients derived for up 
to the 4-th harmonic. The results were a bit disappointing 
in that all types of formulae yielded rather large standard 
deviations and of similar order V w 0.1 mag. In Table 4 we 
list selection of the best formulae for each length category. It 
must be kept in mind, that according to the F test, none of 
these formulae performed significantly better than the rest, 
at 0.95 level. Note that the phase term present in formulae 
F2 and F3 does not vary by more than 0.1 mag. For all these 
reasons we recommend the use of the simplest Fl formula: 



M v = -0.259 - 1.338 log P - 0.726Ai 



(12) 



The corresponding My — < V > relations shown in Fig. 
6 reveal that the predicted My vary across fewer than half 
of the observed range. It is rather disturbing to note that 
all stars except for the 3 most outlying ones form a broad 
horizontal clump consistent with no correlation of My with 
V. It is hard to imagine a more vivid demonstration that 
the observed luminosity and shape of the RRc light curve 
do not follow as tight a relation as the theoretical one (see 
point iv in Sect. 3.1.1). 

3.2.2 RRc as standard candles in LMC 

The question whether the shape of the RRc light curve is 
or is not correlated to the luminosity in principle does not 
exclude their use as standard candles. More problems in this 
respect stem from the large intrinsic scatter of their mag- 
nitudes, of order 0.1 mag. For LMC Zebruri et al. (2001) 
list OGLE observations of 450 RRc stars from the fields 
SC2, SC3, SC4 and SC5. As these observations were only 
a byproduct of a project designed for other purposes, the 
number of V filter observations, their exposure length and 
distribution are not optimal for measuring such fine effects 
as the third harmonic in the small amplitude, near sinusoidal 
light curves of RRc. For this practical reason we could not 
test Simon & Clement (1993) formula (Eq. 11) as none of the 
OGLE light curves yielded <^3i with the required accuracy 
of 0.2 radians. 

Application of our formula (Eq. 12) for these data was 
straightforward, as finding the amplitude of the OGLE V 
light curves posed no difficulty at all. From the total sample 
of 450 RRc stars we selected only 57 variables with suf- 
ficient amplitude A\ > 0.1 mag and errors not exceeding 
0.010 mag. In Fig. 7 we plot My computed from Eq. (12) 
against the observed average magnitude. In this respect we 
were reassured, in a perverse way, that the LMC RRc stars 
behave quite similarly to our stars from uj Cen: they both 
reveal little correlation of the predicted My with its actual 
value. Such a situation calls for a repeated analysis of the 
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Figure 7. The My- < V > relation for the 57 OGLE RRc 
variables from the LMC. The solid line has a slope of unity. 
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Figure 8. The A\ — P relation for the 57 OGLE RRc variables 
from the LMC (open circles) and 55 CASE RRc variables from u> 
Cen. 



type performed by KW01 for as large a sample of RRc stars 
from different globular clusters as possible, in order to check 
whether the formulae based on Fourier coefficients have any 
prediction value for My. The impression of similar prop- 
erties is further confirmed in the P — Ai diagram (Fig. 8) 
plotted for RRc stars from both u> Cen (full circles) and 
the LMC (open circles). Both samples occupy the same lo- 
cations without showing any systematic differences. In our 
opinion, this particular similarity justifies use of the RRc 
stars as standard candles in these clusters. 



4 DISCUSSION 

The F-band photometry of stars in ui Cen obtained by 
Kaluzny et al. (2003) contains 3 times as many observations 
and with higher precision than available so far. To under- 
stand what could be learned from such data we attempted 
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Table 3. Basic parameters of RRc in u> Cen based on the photometry of Kaluzny et al. (1997) and Kaluzny ct al. (2003). 
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here to redo from scratch analysis of the RR Lyr stars in 
this cluster. The results combined with the new metallicity 
determinations of Rey et al. (2000) and distance modulus 
to the cluster derived by Kaluzny et al. (2002), yield for RR 
Lyr stars a new Mv — [Fe/H] relation independent of previ- 
ous ones. For evolved stars, with P > 0.7 days the relation 
becomes nonlinear. Rejecting them, we obtain a nearly lin- 



ear relation with inclination consistent with previous results. 
Applying this relation to available observations of RR Lyr 
in the Large Magellanic Cloud we obtain its distance mod- 
ulus (j,lmc = 18.56 ± 0.14 mag consistent with the "long" 
distance scale. 

Using the same V-band photometry we derive and cali- 
brate the formulae connecting the absolute magnitude of RR 
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Table 4. Formulae for My of RRc stars with their T> parameters and the resulting distance moduli to ui Cen with its standard deviations 



No. 


Equation 






X> 






1 


My = -0.259 


- 1.338 log P + 0.726^1 




0.0998 


14.0883 


0.0140 


2 


M v = -0.318 


- 1.330 log P + 4.620A 3 - 


h 0.042(^21 


0.0987 


14.0875 


0.0138 


3 


My = -0.653 


- 2.025 log P - 6.981A 2 - 


h 14.639A 3 + 0.048</-2i 


0.0954 


14.0892 


0.0131 



Lyr stars of Bailey type ab with the Fourier parameters de- 
scribing their light curves. This approach, closely following 
previous works, proved only partially successful in that the 
residuals from fit remain large compared to the statistical 
errors of observations and in that the new formulae show lit- 
tle relation to the ones derived from older data. Employing 
our best fit formula calibrated with our distance modulus of 
ui Cen to the OGLE photometry of 13 RR Lyr stars from 
the LMC enabled us to obtain its distance modulus equal to 
Lilmc = 18.51 ± 0.07 r ± 0.04 s mag. 

The agreement with the value obtained from the My — 
[Fe/H] relation is very good, which on one hand is hardly 
surprising as both data sets overlap, yet on the other hand 
demonstrates the formal validity of our methods. Summa- 
rizing, we estimate the LMC distance modulus equal to 
Hlmc = 18.52 ± 0.06 r ± 0.04 s mag where the systematic 
error is determined by the error in the distance modulus of 
ui Cen (Kaluzny et al. 2002). The distance to the LMC is 
thus equal to 50.6 ± 1.6 kpc. 

Our results for the RRc stars are less encouraging. De- 
spite nearly sinusoidal light curves, our photometry proved 
to be accurate enough for analyzing the harmonics, however 
any fitted formulae employing Fourier coefficients yielded 
Mv only loosely correlated with < V >. On the theoreti- 
cal side this indicates dependence of < V > on some ad- 
ditional factor yet to be identified. On the practical side 
any estimated My suffer from large random errors. A com- 
parison with the RRc stars from the LMC observed by the 
OGLE team (Zebruri et al. 2001) proved difficult, because 
apart from the intrinsic errors of our formulae, the poorer 
sampling of data, obtained for other purposes, prevented an 
accurate evaluation of the harmonics. 

An excellent review of the recent distance determina- 
tions to the LMC was given by Walker (1999). Here we 
point out only the discrepancy between so-called "short" 
and "long" distance scales. The former, based on the sta- 
tistical parallax and the red clump methods, gives the dis- 
tance modulus to the LMC as approximately 18.2-18.3 mag 
(Gould & Popowski 1998, Udalski 2000a, 2000b). This is 0.2- 
0.3 mag smaller than the value of about 18.50 mag resulting 
from the cepheid period-luminosity relation and theoretical 
models of the horizontal branch or globular clusters main 
sequence fitting. 

However, recent improvements of the red clump method 
based on the infrared photometry seem to indicate that this 
method also provides the distance modulus of the LMC 
at around 18.50 mag. For example, three papers based on 
the Tf-band photometry were recently published giving the 
distance modulus to the LMC as equal to 18.49 ± 0.04, 
18.54 ± 0.10 and 18.501 ± 0.008 r ± 0.045 s mag, respectively 
(Alves et. al 2002, Sarajedini et al. 2002, Pietrzyriski & 
Gieren 2002). 

Also, the newest determination of the cepheid period- 



luminosity relation, based on ~600 stars, indicates that the 
distance to the LMC is around 18.50 mag (Sebo et al. 2002). 

Our estimate is in excellent agreement with these recent 
results proving that our calibration of the absolute magni- 
tudes of RR Lyr stars produces valuable results. 
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